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NONAXISYMMETRIC INCOMPRESSIBLE HYDROSTATIC 


PRESSURE EFFECTS IN RADIAL FACE SEALS 
by Izhak Etsion* 

Lewis Research Center 

SUMMARY 

A flat seal having an angular misallnement is analyzed, taking into account the 
radial variations in seal clearance. An analytical solution for axial force, tilting mo- 
ment, and leakage is presented that covers the whole range from zero to full angular 
misallnement (surfaces in contact). It is shown that nonaxisymmetric hydrostatic 
pressures due to the radial variations in the film thickness have a considerable effect 
on seal stability. When the hi^ pressure is on the outer periphery of the seal, both 
the axial force and the tilting moment are nonrestoring. This causes the seal surfaces 
to wear at the outer diameter. Instability and wear at the inner diameter can occur 
when angular misallnement is combined with radial distortions and the high pressure 
is on the inner periphery. 

The case of high-pressure seals where cavitation is eliminated is discussed, and 
the possibility of dynamic instability is pointed out. 


INTRODUCTION 

Axial forces in excess of those theoretically predicted for alined flat seal faces 
are now well established by many es^eriments (e. g. , refs. 1 to 3). These axial forces 
are crucial in providing the mechanism that separates the seal faces and controls the 
gap between them. 

A considerable amount of fundamental research has been done in the last two 
decades in an attempt to understand the origin of these axial forces and how they relate 
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to the mechanism of seal operation. Many theories have been suggested, but there is 
no generally accepted mechanism. In a recent report (ref. 4), the various models and 
theories are listed. It seems that most of them relate the separating axial force to 
hydrodynamic effects. 

To generate hydrodynamic pressures with relative circumferential motion, the 
geometry of the seal faces must deviate from being absolutely alined or perfectly flat. 
This deviation can be in the form of angular misalinement (ref. 5), microasperities 
(ref. 6), and various forms of waviness (refs. 2, 7, and 8). It was postvilated (ref. 9) 
that cavitation must be associated with the production of net forces greater than the 
hydrostatic forces. In the absence of cavitation, a hydrodynamic pressure increase in 
the circumferentially converging portions of the seal clearance is covinterbalanced by 
an equal pressure decrease in the diverging portions, with the result of zero net hydro- 
dsmamic force. Although it explains the existence of hydrodynamic generated forces, 
cavitation is not likely to occur in hi^-pressure seals, hi some experiments (refs. 

10 and 11) where regions void of fluid were observed, the applied pressure was usually 
less than 3 atmospheres. In reference 11, where temperatures in the seal clearance 
were measured, it was found that boiling rather than cavitation was the reason for the 
appearance of gas in the seal interface. 

Considering all these facts, it is not likely that hydrodjmamic effects are the only 
mechanism that is responsible for maintaining a lubricated gap in face seals. Some 
investigators (refs. 1, 12, and 13) studied the effects on seal behavior of hydrostatic 
pressures caused by radial variations in seal clearance. Such variations, which can 
be caused by thermal and mechanical distortions, generate a net force that is different 
from the hydrostatic one that is related to flat, alined seal faces. 

When the clearance decreases in the direction of flow, the pressures are greater 
than those for parallel surfaces. Also positive stiffness and hence a stable fluid film 
are related to such configurations. When the clearance increases in the direction of 
flow the result is a decrease in pressure, negative stiffness, and a statically unstable 
fluid film. 

In all the literature dealing with hydrostatic effects, with the exception of refer- 
ence 14, the deviation from the flat, alined geometry is always axisymmetric. The 
inevitable result is axisymmetric pressure distribution in the seal clearance. Some 
experimenters report surface rubbing in a form that is most unlikely to occur when 
pressure symmetry prevails. In reference 1 it is shown that the seal surfaces rub on 
either the outer or the inner boundaries, depending on the direction of flow. The wear 
is always such that film convergence in the direction of flow is produced. Obviously, 
such wear could be a result of an angular misalinement combined with nonrestoring 
moments. 
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In reference 14 a misalined hydrostatic seal with an orifice bypass is analyzed. 
For the case of zero orifice flow, which is no more than a conventional face seal, it is 
found that the seal always produces restoring tilt moments. This result contradicts 
the experimental evidence of surface wear (ref. 1), which suggests nonrestoring tilt 
moments. 

Considering the misalined face seal of figure 1, it is clear that, viiile the clear- 
ance converges in the radial direction over one-half of the circumference, over the 
other half it diverges. Hence, when the flow is outward, the hydrostatic pressure dis- 
tribution around the seal tends to reduce the tilt and restore alinement. On ihe other 
hand, when the flow is inward, the result is a nonrestoring moment, which tends to in- 
crease misalinement. 

Nonaxisymmetric hydrostatic pressure distribution may affect not only the tilting 
moment, but also the separating force and the leakage across the seaL Such effects 
have been overlooked in other investigations. The objective of this work is to study 
these effects, aiming at a better imderstandlng of radial face seal operation. 
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SYMBOLS 

seal clearance along centerline 
constants of integration 
axial force 
hydrod 5 mamic force 

axial force contributed by nonaxisymmetric hydrostatic pressure 

nondimensional force, F/ 7 r(pj - p^)r^ 

nondimensional film thickness, h/C 
film thickness 
given by eq. (15) 
tilting moment 

O 

nondimensional tilting moment, M/ir(p. - Pq)i*q 

pressure 

leakage 


3 



Q nondimensional leakage, Q / 


c- /Pi - Po 


\ r. - 


+ Tj) 


R nondimensional radius, r/r^ 

r radial coordinate 

/3 angle of distortion 

y angle of tilt 

e tilt parameter, yr^/C 

0 angular coordinate 

fx viscosity 

CO rotational angular velocity 

Subscripts: 

1 at inner radius 

m at midradius 

o at outer radius 


ANALYSIS 

Figure 1 describes the geometry of the misalined seal. In most applications the 
radius ratio is very close to unity. Hence, the narrow- seal approximation re- 

sults in the one- dimensional incompressible Re 3 molds equation 

/rh^ = 6^^(or (1) 

9r y drj dO 

The film thickness distribution for the misalined seal is given by 

h = C + yr cos 9 (2) 

where y is the angular tilt and 9 is measured from the point of maximum clearance. 
The boundary conditions for equation (1) are 
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p = p. at r=r. 

P = Po '' = ''0 

Because of its linear nature, equation (1) can be solved separately for the hydrostatic 
and hydrod 3 ntiamic pressure components. Hence, the hydrostatic component is obtained 
from the solution of 

(rh^ = 0 (3) 

8r \ 8r/ 

A common assumption for narrow seals is that curvature effects may be neglected. 
When applied to equation (3), this must be carefully examined. In some works tiiis as- 
sumption results in a modified ejqpression for the film thickness having the form 

h = C + vr^ cos 0 
' m 

This is, of course, an unrealistic, warped surface that if substituted into equation (3) 
will result in hydrostatic pressure distribution identical to that in the case of parallel 
seal surfaces. 

Hence, when neglecting curvature effects, equation (3) becomes 

A /h^^^ = 0 (4) 

8r \ 8r/ 

but the film thickness h remains a function of both 6 and r, as given by equation (2). 
Integrating twice with respect to r gives for the pressure 

C,(0) 

P=^— +C (0) 

V.2 2 


For the boundary conditions, the constants of integration become 


and 


Ci=(Pi-Po) 


<W 

v,2 u2 

^o" h 
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Hence, the hydrostatic pressure in the seal clearance is 


/h.\^ h^ - h^ 
P=(Pl-Po)(^] 


( 5 ) 


Using the film thickness expression given by equation (2) yields 


h2-h2 




h2-h2 ^o-^i^o-^^i 


Substituting this into equation (5) yields 


P = (Pi - Pq) — 1 + 

^o-^i 



h. + h 
o 

h +h. 
o 1 


- 1 


+ Pr 


(6) 


The first term in the braces of equation (6) represents the well-known axisymmetric 
solution for parallel surfaces. The second term is the contribution of misalinement. 
This contribution is not axisymmetric and depends on the gap. 


Axial Force 

The axial force, which tends to separate the seal surfaces, is obtained from 


F = 2 


r° r 

-'r. ''0 


pr dr d0 


( 7 ) 


The nondimensional film thickness 


H = ii- = 1 + eR cos 0 
C 


is introduced where 


( 8 ) 
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and the tilt parameter e is given by 


e = 



Equation (5) can be written in the form 


P = (Pi - Pq) 




+ Pr 


( 9 ) 


A direct integration of equation (7) Is very difficult. A perturbation technique (ref. 5) 
could be used where the pressure is expanded in powers of e. However, such a solu- 
tion is restricted to very small values of the tilt parameter e and is inadequate when 
nonrestoring moments exist since this situation may lead to surface contact. In that 
case the solution has to cover the whole range of misalinement through e = 1. 

With the assumption of negligible curvature effect, that is, r « r^^^, equation (7) 
becomes 


= 2r r n p dR d0 

TO. O Jr. Jd 


0 •'R. 


Integrating by parts we have 


rl 1 rl 

j p dR = pR - / 

\ R. \ 


R^ dR 
0R 


From equations (9) and (8) we have 

_ 9E.= 


oti2tj2 

, , '^"i % e cos 9 

(Pi - Pq) 7~ 

0R 


Substituting equation (12) in equation (11) and using the integral 


/ 


eR cos 6 ^ 1 

3 


1-1+^ 


H" 


e cos e \ H 2 h^ 


( 10 ) 


( 11 ) 


( 12 ) 


we have 
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( 13 ) 


where 


Since 


and 


/ p dR = Pq - PjR. + (Pj - Pq)I 

•'R. 


1 , 2H.H 

I = — ^ { i-2-- 1 

e cos 6 \ Hq + H. 


Ho ^ Hi = 


1 + R. = 2R^ 
1 m 


the expression for I becomes (by using eq. (8)) 

R + eR. cos 6 


1 = 


_ m i 


1 + e cos 6 


where 


The axial force thus becomes 


e = eR. 


m 


F = 2r^r. 
m o 


7T(P, 


o - PiR;) + (Pi - P„) I 


dd 


From the journal bearing integrals (ref. 15), 

► 7T 


/ 


de 

1 +e cos 6 


IT 


- e^) 


. 1/2 


and 


(14) 


(15) 


(16) 


(17) 


(18a) 


8 




cos 0 dg 
1 +'i cos d 


Hence, from equations (15) to (18) 



(18b) 


r 


F = 27rr r J 
m o> 


Po - pA + (Pi - Pq) 


Rtv, - 
m 


h. 


/ 2 2 

(1 - ^ Km) 


m 








(19) 


This can be further developed into 


r 


R.. 


^ - V Pi^^m - ^i) (Pi - Po) ( ^m - ^ 








1/2 


- 1 




Noting that 


and that 


1 - 

m m 


R = 1 (1_ R) 
1 2 ^ 


R. 


^m- 


m 



(1 - Ri)^ 


4R 


m 


(20a) 


(20b) 


we finally have for the axial force 


^ = ’^m(''o - ""iXPi + Po> ’"''o(Pi - Po> 




1/2 


- 1 


( 21 ) 


The first term on the right side of equation (21) is the axial force for the case of 
parallel flat surfaces. The second term presents the contribution of the nonaxisym- 
metric hydrostatic pressure distribution due to seal misalinement. Denoting that con- 
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tribution by F , we have 

s 

F = 7rr^(p. - pJF 
s o 1 ^o' s 

where F is a nondlmensional force given by 
s 

? 

»r2(pj - 2 

From equation (22) it is clear that, for any angular misalinement, F^ is positive. 
Hence, the deviational force F^ depends on p. - p^. For p. > p , F is positive 
and thus yields axial forces in excess of the hydrostatic force for alined flat surfaces. 
Since e is given by yr^/C, it can also be seen that a decrease in C increases e 
and hence increases F^. Thus, when p, > p_, the axial force F_ increases when the 
clearance decreases, producing positive stiffness and a statically stable film. On the 
other hand, when p. < p^, the axial force F^ and the axial stiffness become negative 
and thus result in vmstable operation. 




- 1/2 


(22) 


Tilting Moment 

The tilting moment due to the nonaxLsymmetric component of the hydrostatic pres- 
sure is 


M = - 



cos Q dr d0 


The (-) sign is used in order to obtain restoring moments as positive values. Again, 
when the curvature effect is neglected, the moment becomes 


M = -r^r 


r f 

® *^R. *4) 


The same procedure as for the axial force is 
in a form similar to equation (17): 


“ = -2 Vo'Pi - Po> 


27T 

p cos 6 dR 6.6 (23) 

repeated here, giving the tilting moment 



I cos 6 d0 


(24) 
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From reference 15 we have 



cos^.^._d0=JL 

1 + e cos 6 —2 




1/2 


- 1 


(25) 


Substituting equation (15) into equation (24) and using equations (16), (18b), and (25) 
give the tilting moment as 


M= (Pj - PQ)7rrQM (26) 

where the nondimens ional tilting moment M is given by 

As in the case of the axial force, the nondimensional moment M is always positive 
and vanishes only for e = 0, as expected. For p. > p^ the tilting moment as given by 
equation (26) is positive and hence is a restoring moment. For p. < p^ the tilting mo- 
ment is a nonrestoring one, the angular stiffness is negative, and the seal is unstable. 


- (1 - RJ' 

M = 1- 


1 

6 



Leakage 


The leakage across the seal is obtained from 


Q = 2 


r_ hi 

Jo 12m 


^ rde 
9r 


(28) 


If we use the nondimensional terms H and R, substitute equation ( 12) into equa- 
tion (28), noting that - H? = (H^ + Hj)(l - R.)e cos 9 , and neglect the curvature ef- 
fect, equation (28) becomes 


Q 




cos 6 


de 


(29) 


11 


H?Hq = 1 + 2(1 + R.)e cos 0 + 1^(1 + + 2RjJe^ cos^e 


+ 2Rj(l + Rj)€^ cos^e + Rpe'^ cos‘^0 


Again, from reference 15 we have 


and 


I 

I 


_22Ll_d0^X 


I + € COS 0 —3 

e 


1 - 


L k--^) 


1/2 2 


4« 

cos 0 d0 = — 


1 + € COS 0 —4 

c 




1/2 


- 1--^ 


Using equations (18), (25), (31), and (30), we have from equation (29) 

r 


/ R 

C / m 








+ 4 


^i 


r: 


m 



— 

p— . 



1 

1 


< 


)i/a ^ 


[_ 1 

-2 

4- e 


1 



L(i - ^ ) 


After rearranging this becomes 




12k \r^- r, 


where the nondimensional leakage Q is given by 


(30) 


(31a) 


(31b) 



■JJ 


( 32 ) 
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From equation (20b). Q can be written in the form 



From equation (33) it is clear that for the alined seal, where e 0, Q will equal 1, 

as expected. As the misalinement increases, tiie leakage increases too. For very 

narrow seals, where 1 - R. < 0. 1, it is seen from equation (33) that the increase is 

1 2 

most affected by the term (3/2) e; The nondimensional leakage, given in the form 

Q-l+-e^R,. (34) 

2 ^ 

is obtained when the radial variations in the film thickness h (eq. (2)) are neglected. 
Hence, it can be seen that, with regard to leakage, the simpler expression 
h = C cos 6 can be used in the analysis. 

RESULTS AND DISCUSSION 

Values of the nondimensional parameters F_, M, and Q are presented in table I 
and figures 2 to 4 to cover the whole range of tilt parameters from e = 0 to e - 1. 

An interesting result is the axial force at e = 1. In the case of radially inward flow, 
this force adds to any imbalanced closing force and Increases substantially the friction 
between the rubbing surfaces. A seal of 10- centimeter outer diameter, 0. 9 radius 
ratio, and 10-atmosphere seal pressure, will experience at e - 1 an unbalanced axial 
force of 86 newtons and a tilting moment of 430 newtons per centimeter. 

The radius ratio r./ r^ affects very much the axial force and tilting moment. As 
this ratio decreases, the axial force and tilting moment increase quite rapidly. 

As can be seen from table I the differaices between the values of Q calculated by 
equation (33) and those calculated by equation (34) ai*e negligible. Hence, the simpler 
egression (eq. (34)) can be used for any practical radius ratio. 
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The analysis presented in tlie previous section has indicated unstable seal opera- 
tion under radially inward flow when pj < p^^. Sucli instability' is confirmed by tlie ex- 
periments of Denny (ref. 1). The nonrestoring moment in the case of i-adially inward 
flow causes the surfaces to contact and wear at the outer diameter of the seal. How- 
ever, Denny also found that with the flow radially outward the wear was greatest at the 
inner diameter. Such a wear appearance is unlikely for a flat misalined seal but can 
be a result of combined mislincment and coning due to thermal and mechanical distor- 
tions. 

In figure 5, a misalined distorted seal is shown. It is clear from the figure that, 
whenever the distortion angle /I is larger than the tilt angle y, the seal surfaces can 
come into contact at the inner diameter. Moreo\ er, since an unstable film is related 
to diverging clearances in the direction of flow, it can be seen that, when P > y, an 
instability will result from a radially outward flow. 

A word should be added with regard to the effect of the hydrodynamic pressure 
generated due to angular misalinement. Figure 6 shows the resultant hydrodynamic 
force for a misalined seal. While the hydrostatic pressure distribution is sym- 
metric about the line BB and therefore produces moments only about the line AA, this 
is not the case with the hydrodymamic pressure. If the sealed pressure is high enough 
to eliminate any cavitation, the hydrodynamic pressure distribution becomes antisym- 
metric about the line BB. This results in a net zero hydrodynamic axial force but pro- 
duces a couple that lags the hydrostatic tilt moment vector by 90°, Such a system is 
prone to dynamic instability that can damage the seal. 


CONC LUDING REMARKS 

Contrary to the common belief that hydrodynamic effects are mainly responsible 
for proper seal operation, it is shown that hydrostatic effects play an important role 
in seal performance. The hydrostatic pressure distribution is very sensitive to radial 
variations in film thickness. Hence, these variations cannot be neglected even when 
dealing with radius ratios close to unity. 

The analsdical solution for a flat radial face seal having an angular misalinement 
shows that hydrostatic axial forces and tilting moments can be generated because of 
the misalinement. These forces and moments can be either restoring or nonrestoring, 
depending on the location of the high pressure with respect to the seal boimdaries. 
When the high pressure is on the inside, the seal is stable. When the hi^ pressure is 
on the outside, static instability occurs. 
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The stability condition can be reversed when coning, due to thermal and mechani- 
cal distortions, is combined with angular misalinement. Finally, hydrodsmamic pres- 
sures due to seal misalinement can well be the origin of dynamic instability. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 30, 1976, 

505-04. 
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Ratio of 
inner to 
outer 
radiiis, 

O.BO 


0.8S 


0.90 


0.92 


TABLE I. - MISALINED-SEAL PERFORMANCE PARAMETERS 


Tilt 

parameter, 

€ 

Nondi mensional 
force, 

Nondimensional 
tilting moment, 
M 

Non- 

dimensional 

leakage, 

Q 

Simplified 

nondimensional 

leakage, 

2 * 

O.IOOOE 00 

0.f^l4«6-04 

0*81486-03 

1,012 

1,012 

0.2000 

0,33216-03 

0.15506-02 

1,048 

1,043 

0.3000 

0,77146-03 

0.25716-02 

1,108 

1,108 

0,4000 

0,14376-02 

0,35936-02 

1,193 

1,102 

0,5000 

0,23956-02 

0.47916-02 

1,301 

1,300 

O.GOOO 

0*37626-02 

0.52716-02 

1,434 

1,432 

0.7000 

0.57536-02 

0.82196-02 

1,591 

1,583 

0.0000 

0.R0196-02 

0.11026-01 

1,771 

1,763 

0*no00 

0.14106-01 

0.15576-01 

1,976 

1,972 

1,0000 

0.25RR6-01 

0*25886-01 

2*206 

2,200 

o.inooE 00 

0.40436-04 

0.48436-03 

1,013 

1,013 

0,2000 

0.19756-03 

0,98806-03 

1,051 

1.051 

0.3000 

0.45996-03 

0.15336-02 

1,115 

1,115 

0,4000 

0,05946-03 

0.21486-02 

1,204 

1,204 

0,5000 

0. 14396-02 

0. 28776-02 

1,319 

1,319 

O.fiOOO 

0,22746-02 

0,37906-02 

1,450 

1,459 

0.7000 

0,35136-02 

n,50186-n? 

1,525 

1,525 

0,0000 

0.54756-02 

0,58456-02 

1,818 

1,81C 

0,9000 

0.90556-02 

0.10056-01 

2,0’;5 

2,033 

i.onoo 

0,10356-01 

0.18356-01 

2,278 

2,275 

O.IOOOP 00 

0,22716-04 

0.22716-03 

' 1 

1,014 

1,013 

0.2000 

0.92776-04 

0.45386-03 

1,054 

1.054 

0,3000 

0,21536-03 

0.72116-03 

1,122 

1,1.21 

0.4000 

0.40556-03 

0.10146-02 

1,215 

1,216 

0.5000 

0,58196-03 

0.13546-02 

1,338 

1,337 

0,9000 

0.10856-02 

0,18096-02 

1,486 

1,486 

0.7000 

0,15956-02 

0,24216-02 

1.662 

1,551 

0,0000 

0.25936-02 

0.33576-02 

1,865 

1,RG4 

0,9000 

0.45416-02 

0,51556-02 

2,0 95 

2,093 

1.0000 

0.11016-01 

0.11016-01 

2,351 

2,350 

O.lOOOli 00 

0.14856-04 

0,14856-03 

1,014 

1,014 

0.2000 

0.50556-04 

0.30336-03 

1,055 

1,055 

0.5000 

0.14156-03 

0.47196-03 

1,124 

1,124 

0,4000 

0.25576-03 

0,55436-03 

1,221 

1.221 

0.5000 

0.44776-03 

0.89546-03 

1,345 

1,345 

0.6000 

0,71456-03 

0,11916-02 

1,497 

1,497 

0.7000 

0.11216-02 

0*15026-02 

1,577 

1,575 

O.SOOO 

0,17956-02 

0,22456-02 

1,884 

1,883 

0,9000 

0.31556-02 

0,35056-02 

2,118 

2,1.13 

1.0000 

0.82286-02 

0,82286-02 

2.381 

2,3^0 


17 



TABLE I. - Concluded. 


parameter, 


Nondimensional 

force, 


Nondimensional 
tilting moment, 


dimensional 

leakage, 

Q 


Simplified 
no ndi mens ional 
leakage, 


O.IOOOE 00 
0,7.0 00 
0,3000 
0,4000 
0,5000 
0.6000 
0.7000 
0,8000 
0,0000 

1,0000 


0,n528E-05 
0,34866-04 
0,81426-04 
0.15306-03 
0. 25836-03 
0.41356-03 
0.65106-03 
0.10546-02 
0.18016-02 
0,56046-02 


0,85286-04 

0.17436-03 

0.27146-03 

0,38256-03 

0,51666-03 

0.68026-03 

0,03126-03 

0.13176-02 

0,21016-02 

0.56046-02 


O.IOOOE 00 
0,2000 
0,3000 
0.4000 
0.5000 
0,6000 
0,7000 
0.8000 
0,0000 

1,0000 


0,38696-05 

0,15826-04 

0.36996-04 

0.69606-04 

0,11776-03 

0.18906-03 

0,29956-03 

0,48876-03 

0,89766-03 

0.32206-02 


0.38696-04 
0,79126-04 
0,12336-03 
0.17406-03 
0,23546-03 
0. 31516-03 
0,42796-03 
0.61096-03 
0,99736-03 
0,32206-02 


O.IOOOE 00 
0.2000 
0. 3000 
0,4000 
0.5000 
0.6000 
0.7000 
0,8000 
0.9000 

1,0000 


0.98726-06 

0,40396-05 

0,94516-05 

0.17816-04 

0.30186-04 

0.48616-04 

0.77426-04 

0,12766-03 

0,24056-03 

0.12186-02 


0.98726-05 
0. 20206-04 
0.31506-04 
0,44516-04 
0.60366-04 
0,81026-04 
0.11066-03 
0,15956-03 
0.26726-03 
0.12186-02 


0.10006 00 
0 ,2000 
0.3000 
0.4000 
0.5000 
0.6000 
0. 7000 
0,8000 
0,9000 

1,0000 


0,24936-06 

0,10206-05 

0.23886-05 

0.45036-05 

0.76396-05 

0,12336-04 

0.19686-04 

0.32606-04 

0.62346-04 

0.45066-03 


0,24936-05 

0.51026-05 

0.79616-05 

0.11266-04 

0,15286-04 

0.20546-04 

0,28126-04 

0.40766-04 

0.69276-04 

0.45066-03 
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